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I review recent developments in Heavy Quark Effective Theory (HQET) that lead 
to an almost model-independent determination of the \ V c b\ element of the Cabibbo- 
Kobayashi-Maskawa matrix from exclusive semileptonic B — > Z)M decays. In par- 
ticular, I compare the theoretical uncertainties in the B ^ D*£ u and the B D £u 



Abstract 



decay modes. I discuss the applications of QCD sum rules within HQET to semilep- 
tonic heavy meson decays and give predictions for the form factors measurable in 
B -> D&lv decays. 



★ Invited talk at the Advanced Study Conference on Heavy Flavours; September 3-7, 1993 at Pavia, Italy. 



1. Introduction 



In the absence of direct observations of new physics, heavy hadron decays may provide 
the first clues to physics beyond the Standard Model (SM). They probe the flavor sector, 
which contains the majority of the parameters of the SM. These parameters have to satisfy 
certain relations provided by, e.g., the unitarity of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix, and the SM prescription of CP violation. Measuring the parameters of the SM more 
accurately, one hopes to find inconsistencies among them, which would give hints to new physics. 
Even if the numerical values of these parameters are consistent with each other, they may 
provide insights to new physics that can yield relations among them. (For example, schemes 
for quark mass matrices can be tested, that may teach us about horizontal symmetries or grand 
unification.) Certain rare decays are particularly sensitive probes of various extensions of the 
SM. However, the theoretical predictions for most of these measurements are contaminated 
by large hadronic uncertainties. Reducing these uncertainties would provide better chances of 
discovering physics beyond the SM, and make the bounds on such new physics more restrictive. 
Recent developments in this direction constitute the subject of this talk. 

In hadrons composed of one heavy quark and a number of light degrees of freedom (gluons 
and light quarks), the energy scale of strong interactions is small compared to the heavy quark 
mass. The heavy quark acts effectively as a static source of color, resulting in new symmetries 
of QCD [1-7]: the interaction between the heavy quark and the surrounding light degrees of 
freedom become independent of the mass and spin of the heavy quark. In the tjiq ^> Aqcd 
limit the velocity of the heavy quark is conserved with respect to soft processes [1,7], and the 
complexity of hadronic dynamics results from the strong interactions among the light degrees 
of freedom only. From the phenomenological point of view, such a symmetry (even if broken 
as the 'heavy' quarks are, after all, not infinitely heavy) is an extremely useful tool, as it 
provides exact predictions in the symmetry limit, reducing hadronic uncertainties and model- 
dependence, which can only enter in corrections, suppressed by powers of 1/mg. 

The heavy quark effective theory (HQET) [5-11,13] provides a convenient framework to 
analyze heavy hadron decays. It allows for a systematic expansion of hadronic quantities in 
powers of Aqcd/^q i n such a way that the coefficients are heavy quark spin- and mass- 
independent universal functions of the kinematic variable y = v ■ v', where v and v' denote 
the velocities of the initial and final hadrons. These universal functions originate from long 
distance QCD, so they can only be investigated using nonperturbative methods. Such a method 
is provided by QCD sum rules [14-16], which have been widely used recently to calculate 
hadronic matrix elements in HQET [17-24]. 

In this talk I shall focus on decays of heavy mesons (rather than baryons). Experimentally 
they are easier to measure, resulting in more phenomenological applications (e.g., extraction 
of \V c f)\), and most calculations of HQET form factors using either QCD sum rules or other 



models have been carried out for these decays. Moreover, in heavy meson decays a comparison 
with lattice calculations is also possible [27]. 

In Section 2, I review aspects of HQET that are relevant to the discussion in Section 3 
of the almost model-independent determination of \V c i,\ from exclusive semileptonic B meson 
decays. In Section 4, I discuss applications of QCD sum rules for heavy meson decays. Besides 
the universal functions, predictions for heavy meson decay constants and semileptonic form 
factors are also presented. Finally, I summarize and outline some directions of ongoing (and 
future) developments. 



2. HQET 

The construction of HQET starts with removing the mass-dependent piece of the momen- 
tum operator by introducing a field hqiy , x), which annihilates a heavy quark with velocity v [7], 

h Q (v,x)=e im ^ x P + (v)Q(x), (1) 

where P+(v) = ^(1 + t/>) is an on-shell projection operator onto the heavy quark (rather than 
antiquark) components of the spinor, and Q(x) denotes the conventional quark field in QCD. If 
P^ is the total momentum of the heavy quark, the new field Hq carries the residual momentum 
= P^ — itlq ~ V(Aqcd), which does not grow with the heavy quark mass. 

In the limit rriQ ^> AqcDj the effective Lagrangian for the strong interactions of the heavy 
quark is [5-8] 



1 



£hqet = h v iv-Dh v + O kin + C mag (/z) O mag + V(l/m Q ), (2) 

IrriQ L J 

where = — ig s T a Aa is the gauge-covariant derivative. The leading term respects both 
the spin and flavor symmetries. The operators appearing at order 1 / tuq are 

O kin = h v {iDfh v , O mag = | h v a^G^K . (3) 

Here is the gluon field strength tensor defined by [iD^,iD v ] = %g s G liV . In the hadron's 
rest frame, 0^ n describes the kinetic energy resulting from the residual motion of the heavy 
quark, whereas O mag corresponds to the chromomagnetic coupling of the heavy quark spin to 
the gluon field. While Okin violates only the heavy quark flavor symmetry, O mag violates the 
spin symmetry as well. Due to reparameterization invariance [28] Okin is not renormalized to all 
orders in perturbation theory, while O mag (/x) is a renormalization factor for O mag . The heavy 
quark symmetries are also manifest in the Feynman rules of the effective theory: the propagator 
of a heavy quark becomes independent of its mass (flavor symmetry) and no gamma matrix 
appears in the coupling of a heavy quark to the gluon field (spin symmetry). 



f The interested reader can find the HQET formalism for heavy baryon decays in Ref. [25]; the QCD sum 
rules determination of the relevant universal function in [26] . 
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Figure 1: Feynman rules in "full" QCD (a) and in HQET (b). Double lines denote heavy 
quark propagators in HQET (from [29]). 

Any operator of the full theory that contains one or more heavy quark fields can be matched 
onto a short distance expansion in terms of operators of the effective theory. In particular, the 
expansion of the heavy quark current Q' T Q that gives rise to B — > decays is matched 
onto 



2ttiq 2rriQi 



Oj + V(l/m|) . 



(4) 



The operators {Ji} form a complete set of local dimension three current operators with the 
same quantum numbers as the current in the full theory. There are three such operators 
Ji — hQ/TiliQ, with Ti = {7^, yt 1 , v'^} for the vector current, and T; = {7^75, f M 75, 1/^75} 
for the axial current. (In the leading logarithmic approximation only T = 7^ (75) contributes. 
Radiative corrections induce the other operators.) Similarly, {Oj} denote a complete set of 
local dimension four operators. Since there are fourteen independent such operators, we do 
not display them explicitly here. These effective current operators have non-zero anomalous 
dimensions. The coefficients Ci(fi) and Bj('\/j,) ensure that the final result for any physical 
quantity is independent of the renormalization procedure. At present, the expansion of the 
effective Lagrangian and the weak currents are known in perturbation theory up to order 
a s /rriQ and 1/rriQ [13]. 

2.1. Leading order 

Matrix elements in HQET are conveniently calculated in the compact trace formalism [9,10], 
where a heavy meson is represented by its spin wave function 



M(v) = ^Tf^ 




-75 ;Jt 



cr 



; J p = 1-, 



(5) 



which has the correct transformation properties under Lorentz boosts and heavy quark spin 
rotations. When the external weak current changes v — > v' (and maybe Q — > Q'), the light 



degrees of freedom have to rearrange themselves, which yields a form factor suppression. Due 
to heavy quark symmetry, this suppression factor cannot depend on the spin and the mass of 
the heavy quark, neither on the Dirac structure of the current. Lorentz and parity invariance, 
and the properties of JA(v) imply that dependence only on y = v-v' and on the renormalization 
scale ji is allowed. Hence, a single universal, i.e. only y and fi dependent, function £(?/,//) is 
sufficient to parameterize all semileptonic M(v) — > M'(v')£u decays, where M and M' are 
pseudoscalar or vector mesons containing a single heavy quark: 

(M'(v')\h' v ,rh v \M(v)) = -Z(y,ti) Ti{M\v')YM{v)}. (6) 

Vector current conservation implies that when the heavy meson in the final state is at rest 
in the rest frame of the decaying heavy meson, this so-called Isgur-Wise function satisfies 
£(1) = 1. The predictions of HQS are most restrictive at this special kinematic point ("zero 
recoil": y — 1), allowing model-independent predictions, unaffected by hadronic uncertainties. 

Thus, at leading order in the heavy quark expansion, matrix elements factorize into a 
kinematic part that depends on the mass and the spin-parity of the mesons, and a reduced 
matrix element that describes the light degrees of freedom. This is a remarkable simplification, 
as a-priori six independent form factors describe the semileptonic B — > transitions. Since 
the b and c quarks are not much heavier than Aqcd, an analysis of the I/ttlq corrections is 
important for most phenomenological applications. 

2.2. 1/rriQ CORRECTIONS 

At subleading order, matrix elements receive contributions from the higher dimension op- 
erators in the effective Lagrangian (2) and in the effective current (4). The idea is to leave the 
heavy quark propagator identical to its leading order expression and account for the correction 
terms in the Lagrangian as insertions of operators. To parameterize their matrix elements 
we need three new universal functions Xi(y) = 1)2,3). Vector current conservation implies 
= X3(l) — (t ms is known as Luke's theorem [11]). 

Matrix elements of the I/tuq corrections in the effective current (4) are parameterized in 
terms of another three universal form factors, usually denoted by and £3(1/). Im- 

posing the equation of motion, i(v-D)fiQ = 0, on the matrix element yields two constraints [11]. 
Thus only one of these three functions, say £3(1/), is independent. 



function 




xi(y) 


X2(y) 


xs(y) 




normalization 




xi(i) = o 


no 


xs(i) = o 


no 


broken symmetries 


no 


flavor 


spin, flavor 


spin, flavor 


spin, flavor 



Table 1: Properties of the universal functions of HQET. 



So already at order 1/mg one encounters a set of four universal functions £3(2/) and 
Xi(y) = I5 2, 3) in addition to the Isgur-Wise function, as well as a parameter A = uim —mQ 
that describes the mass difference between a heavy meson and the heavy quark that it con- 
tains [11,30]. This parameter sets the scale of the 1/mg expansion; in fact, the real expansion 
parameter is A/2mg. The universal form factors are real due to T invariance of the strong 
interaction. Knowledge of these functions would teach us about confinement and enhance the 
phenomenological applications of the heavy quark expansion. 



3. Model Independent Determination of \V c b\ 

The magnitude of the CKM matrix element can best be determined from an extrapo- 
lation of the semileptonic B decay rate to zero recoil, making use of the known normalization 
of the Isgur-Wise function at that point [31]. The advantage of this method over previous de- 
terminations of \V c t,\ is that in the framework of HQET a clear separation between the model- 
independent and model-dependent ingredients of the analysis is possible, due to a systematic 
expansion in the small parameters A/2m Ci fr . 

The B — > D^iv differential decay rate near zero recoil is given by, 



dT(B^D*£v) G 2 F 3 . . 2 _ |2 . 2 iU/2 , 2 

hm — ^ " = ^3 m D*(^B - m D *) 2 \V cb \ 2 (y 2 - I) 1 ' 2 V * 2 , 



lim ^ B -> Dt *) =§E- 3 ml(m B + m D ? \V cb \ 2 {y 2 - l) 3 ^ 2 . 
2/->l ay 4s7r J 



(7) 



On the right hand sides of these relations the Fermi constant and the meson masses are well 
known quantities, powers of (y 2 — 1) arise from phase space, 7/*) is defined to include all hadronic 
uncertainties at y — 1, and we want to extract \V cb \. The kinematic variable y is related to the 



conventional q 2 via 



, m B + m2 DM ~ g 2 , fi x 
y = v • v = — . (8) 

2m B m D{ ,) 



In this notation maximal q 2 corresponds to y — 1, while q 2 = corresponds to maximal y, 
which is about 1.5 and 1.6 in B — > D* and B — > D decays respectively. The problem is 
that a-priori we know nothing about r]^*\ except that it should be of order one. The power 
of HQS is that it gives the model-independent prediction in the infinite quark mass limit: 
lim mQ ^oo = £(1) = 1. This allows us to write 

77W = 1 + 6$ + 5 { */ mQ + 8 { *] ml + higher order . (9) 

We shall discuss each of the correction terms in the sequel. The calculation of the perturbative 
QCD corrections must include the full order a s terms (not just the leading logarithms), because 



ln(m&/m c ) ~ 1.2 is not a big number. We emphasize that these corrections do not introduce 
uncertainty into the analysis. They are given by 5 as = 0.05 and 5* s = —0.01 [32]. While 
the B — > D*lv decay rate is protected against 1/mg corrections at zero recoil due to Luke's 
theorem [11], i.e. ^[j mQ = 0, the B D£u decay is not, due to its helicity suppression [12,31]. 
In this latter case the 1/mg corrections are given by 

w(£ + £)te^)V*<i>]. m 

Clearly, the form factor £3(2/) is very important for the determination of \V c b\ from B ^ D iv 
decays. For example, if £3(1) were around —1 then this I/ttlq correction would be about 15%, 
while if it were around 0.5 then the 1/mn correction would vanish. In the next section we shall 
see that QCD sum rules predict £3(1) = 0.6±0.2 [24], which implies that the 1/rriQ correction to 
the B — > D tv decay rate at zero recoil is not more than 3%. This is a result of two suppression 
factors (beyond A/2mg): the Voloshin-Shifman factor [(mg — mjj)/(mB + ^z>)] 2 — 0.23 [4], 
and an accidental suppression factor of [1 — 2^(1)] ~ 0.2. The 1/rriQ corrections are expected 
to be 3-4% on dimensional grounds, and the detailed analysis of Ref. [33] supports that these 
corrections are not larger than the above estimate. (This statement, however, is somewhat 
model-dependent, which leaves room for arguments that these corrections might be larger.) 
The higher order corrections are certainly negligible, e.g. the characteristic size of the second 
order QCD corrections is [a s (m c )/7r] 2 < 1%. Thus the 1/rriQ correction to the B — > Dtv 
decay rate at zero recoil is not more than the expected l/m^ corrections. This suggests that 
the theoretical uncertainty in the determination of |V^,| from B — > D transition is comparable 
to that in B — > D*, even though the latter appears only at order I/tuq. Of course, the 
experimental measurement of B — > D tv near zero recoil is more difficult due to extra power 
of (y — 1) helicity suppression in Eq. (7). The reward of such a measurement, however, would 
be an independent determination of |V^,| with surprisingly small theoretical uncertainties. 

A different kind of uncertainty enters into the analysis because phase space vanishes at 
y — 1. Therefore, an extrapolation of the measured spectrum to zero recoil is needed to obtain 
the numerical value of |V^,|. At present this gives rise to both a theoretical and an experimental 
error. The former is due to the fact that the precise shape of the Isgur-Wise function is not 
known, while the latter is dominated by statistical error. However, in the not-so-distant future 
this theoretical uncertainty will almost disappear, since in an asymmetric B factory the zero 
recoil limit does not correspond to the D* meson being at rest in the laboratory frame. Then 
the pion in the subsequent D* — > Dn decay is boosted, while it is almost at rest for ARGUS and 
CLEO. In addition, the (y 2 — l) 1 / 2 phase space suppression is a very mild one: the statistical 
error of measuring the rate at y — 1.05 is less than a factor of two higher than that at the 
endpoint y m&x ~ 1.5. At present any information is important on the shape of the Isgur-Wise 
function, in particular on its slope at y — 1, to make the extraction of |V C &| more reliable. 
Global quark-hadron duality (that the sum of probabilities to decay into hadrons equals to the 



probability of free quark transition when rriQ — > oo) gives the Bjorken [10] and Voloshin [34] 
sum rules: 

\<P 2 ^1, (11) 

where p 2 is minus the slope of the Isgur-Wise function at y — 1. For a discussion of the bounds 
derived by de Rafael and Taron [35] we refer to [36]. 

In summary, we emphasize that due to heavy quark symmetry the theoretical prediction 
for the B — > D^iv decay rate near zero recoil is very precise: 

V* =1 + ^+^=0.99 ±0.05, 
77 = i + 5 as +5 1/rnQ +5 1/ml = 1.05 + 0.08. 

Experimentally only the B — > D*£ V spectrum has been measured (for B — > £) £ P only the 
total decay rate is known). The following figure shows the most recent data from the CLEO 
collaboration [37]. 




Figure 2: The measured dT(B — ► D*lv)jdy distribution. The curves represent fits with 
different functional forms for the Isgur-Wise function. 



They obtain from a fit to this distribution 

\V ch \ =0.037 + 0.005 + 0.004, 
p 2 =1.0 + 0.4 + 0.2. 

As the theoretical uncertainty in the determination of |V^,| is significantly smaller than the 
roughly 15% systematic plus statistical error at present - we should look forward to new data 
from CLEO and a future B factory. 



(13) 



4. QCD Sum Rules 



We have seen, so far, two important examples when more information is needed about the 
universal functions describing heavy meson decays than what is provided by vector current 
conservation (see: Table 1). The shape of the Isgur-Wise function is important for the extrac- 
tion of \V c b\, and knowledge of £3(1) is essential for the determination of \V c b\ from the decay 
B — > D tv '. These universal functions originate from the long distance confining interactions, so 
they can only be investigated using nonperturbative approaches to QCD. While short distance 
QCD is well understood, there is still no accurate quantitative framework for dealing with the 
long distance, strong interaction regime. Lattice gauge theory is the only method known that 
comes from first principles and has the hope of achieving arbitrary accuracy. This is hindered, 
however, at present, by the technical limitations of this approach. A less fundamental method 
is provided by QCD sum rules [14-16], which is particularly suited for the calculation of HQET 
form factors. It has the advantage over other models that it does not rely on ad hoc assump- 
tions, the errors can be estimated, and there is a consistency check on its assumptions within 
the method. 

In the rest of this section we first describe the QCD sum rule method through the example 
of calculating the heavy meson decay constant in HQET, then give the results for the universal 
functions, and finally discuss the predictions of these results for the form factors in semileptonic 
heavy meson decays. 

4.1. Simplest example: decay constant 

The following analysis of the two-point function is particularly important as it determines 
besides the heavy meson decay constant, the parameter A as well. The analogue of the decay 
constant Jm in HQET is defined by 



In leading order fMy/ m M — The idea of QCD sum rules is to calculate the current- 

current correlator in two different ways. In HQET one defines [17] 



Here Jm is an interpolating current for the ground state heavy mesons Jm = h v Ym Q (Tp — —75 
and IY = Ifi — for the pseudoscalar and vector mesons respectively), k = P — rnqv is the 
"residual" off-shell momentum [7], and uj = 2v -k. It is convenient to factor out the Lorentz 



(0\qrh v \M(v)) = ^F(n)Tr{FM(v)}. 



(14) 




(15) 



structure of the two-point function, by defining ir(u) through H(uj) 



±Tr{T M P+T M }Tr(Lo). 



Theoretical side: One side of the sum rule is a diagrammatic calculation of the correlator in 
terms of quark and gluon fields. As we approach resonances from short distances, nonperturba- 
tive effects induce power corrections that violate asymptotic freedom. These effects are taken 
into account in the Wilson operator product expansion (OPE) 

7r the » = C nOn = TTpertM + WlM (16) 
n 

through the so-called "condensates" . These are vacuum expectation values (VEV-s) of gauge 
invariant local quark-gluon operators [14]: 

(0|gg|0), (0\a s G^G^\0), (0\ g s qa^G^q\0) , ... (17) 

The perturbative contributions are depicted in the first line of Fig. 3, while the condensates 
are shown in the second and third lines (the gluon condensate (a s GG) does not contribute). 

-o- 
~~o <?>--- -<^~ 

— ■-~<X > - ^F} C^""" 

Figure 3: Feynman diagrams contributing to the sum rule for meson decay constants in HQET. 
The disjoint lines symbolize the nonvanishing VEV-s of the condensates (from [13]). 

The condensates vanish by definition in the standard perturbation theory. Assuming that they 
can acquire non-zero VEV-s is a simple (but maybe very reasonable) modeling of the QCD 
vacuum. Since the condensates have different dimensions, Eq. (16) forms an expansion in 
inverse powers of uj. The coefficients C n include, by construction, only short distance effects, 
the large distance contributions are accounted for by the vacuum to vacuum matrix elements. 
Strictly speaking, this gives rise to an uncertainty of the method, as we shall assume that 
C n can be calculated perturbatively, while the nonperturbative physics is hidden solely in the 
VEV-s of O n . These parameters are determined from various processes [14,15], 

(qq) = - 0.013 GeV 3 , 
(a s GG) =0.04 GeV 4 , (18) 
(qg a a aP G a Pq) = ml (qq) , m 2 = 0.8 GeV 2 , 

and then their values can be used as inputs in other calculations. Out of these, the gluon 




condensate is the most poorly determined, however, it plays a negligible role in the sum rules 
we shall discuss. This value of a s ^ 0.3 corresponds to the scale /i — 2A ~ 1 GeV, which 
is appropriate for evaluating radiative corrections in the effective theory. Hence we write the 
correlator as 7r thco = 7r por t + 7r con d, and 7r pcr t can further be rewritten as a dispersion integral 



TTpertl^J 



f dv Ppert ^ + subtractions . (19) 
J v — uj — ie 



Phenomenological side: The second calculation of the correlator is in terms of hadrons. H(u) is 
analytic in uj with discontinuities along the positive real axis. The correlator can be expressed 
as a sum over a complete set of physical intermediate states: 

rP he » = y IW0WO)| 2 + subtractions , (20) 
v ' ^ uj x ~ u - ie v ' 

x A 

where the sum is over both discrete and continuum states, and ujx — 2( m X — m Q)- To 
evaluate the right hand side without a detailed knowledge of the spectrum and the matrix 
elements of the excited states, we employ local (in uj) quark-hadron duality. We separate the 
pole corresponding to the ground state meson, and approximate the contribution of higher 
resonances by an integral over the perturbative spectral density above the so-called continuum 
threshold ojq. We find 

^phen.) = J (jj) f ^ Ppert(^) + subtractions (21) 

The QCD sum rules are obtained by matching the theoretical expression with the phe- 
nomenological one: 

- F = l dv Ppert ^ + subtractions + n cond (uj) . (22) 

2A — uo — ie J v — uj — ie 


In order to extract information about the ground state, one has to go to small values of (—ou) to 
enhance the relative weight of the low energy contributions. On the other hand, the theoretical 
calculation is only reliable when (—uj) is large. To achieve a balance between these contradicting 
requirements, an alternative way to probe tt{oj) for small values of {—uj) is by taking derivatives. 
In the limit —uj — > oo and the number of derivatives n — * oo, one is sensitive to the behavior 
of the correlator at scales T = —uj/n. This defines the Borel transformation 

1 a^I uj 71 ( d 



-B™= lim =pr --p . (23) 
T 1 -Lu,n^<x) Tin) \ duj ' 

-uj/n=T 



Employing this operator alters the sum rule (22) in such a way that 



— possible subtraction terms are eliminated; 

— contributions of higher dimensional condensates are factorially suppressed; 

— higher resonance contributions are exponentially damped in the dispersion integral (the 
weight function [v — lu]^ 1 is replaced by exp[— vjT\). 

This yields the final form of the QCD sum rule: 

F 2 ^) e- 2K ' T = J dup peit (u) e~»' T + 7r cond (c). (24) 


To evaluate this sum rule, the continuum threshold ujq, the Borel parameter T, and A have 
to be determined. Notice that by taking the logarithmic derivative of Eq. (24) with respect 
to T _1 we get a sum rule for A. Next, we have to find the continuum threshold c^o such that 
A is stable with respect to variations of the Borel parameter T. Then we use this value of 
A to extract F 2 (fx) from the sum rule (24). The range of c^o and T in which the predictions 
of the sum rule (for A and F 2 (fx) in the present case) are independent of the value of T is 
called the "sum rule window" . This provides a consistency check on the assumption of local 
duality, in the sense that if duality did not hold there were no reason why such a stability 
region existed (it can only result from cancelation of contributions on the theoretical against 
the phenomenological side of the sum rule). For the consistency of the calculation it is necessary 
that both sum rules be stable in the same range of a>o and T. Whether this happens is far from 
trivial and checks, in fact, the consistency of the approximations and the applicability of the 
method. Another assumption of QCD sum rules is that there exists a transition domain where 
the perturbative calculation is still reliable [a s (T) is not too large] and the spectral density is 
sufficiently sensitive to the ground state. This can be assured by requiring that in the sum rule 
window the nonperturbative contributions to the sum rule be less than 30% of the perturbative 
ones, and that the pole contribution accounts for at least 30% of the perturbative part of the 
correlator. 

Several authors calculated the decay constants of heavy mesons [17,18]. The symmetry 
breaking corrections (both from QCD and 1/mg effects) are large, resulting in significant 
deviations from the fs/fD — a/ M^/Mb scaling of the symmetry limit. We quote [17] 

f D = 170±30MeV, fg = 190 ± 50 MeV , (25) 

which compares well with lattice results [27]. The mass parameter A that sets the characteristic 
scale of the 1/mg expansion is significantly higher than Aqcd: A = 0.57 ± 0.07 GeV [17,13]. 



4.2. Universal functions 



For the calculation of the universal form factors, we need to consider three-point correlation 
functions. The calculation, in principle, is very similar to the one outlined above. There are 
certain technical complications related to the applicability of local duality in two variables; 
different choices of the duality region to model the higher resonance states; and the possibility 
of introducing non-local condensates. These have all been discussed in detail in the literature 
[17,19]. The important point is that the physical predictions are obtained by taking the ratio 
of the three-point and the two-point sum rules. This assures that the normalization conditions 
at zero recoil are exactly fulfilled, independent of the value of the Borel parameter T and the 
continuum threshold c^o- These sum rules are probably more accurate than the one for the 
decay constant, since some of the systematic errors drop out from such a ratio, and there is 
no dependence on the less accurately determined A parameter. It is, in fact, very promising 
that the position of the sum rule window is almost identical in the calculations of the decay 
constant, Isgur-Wise function, and the four subleading universal functions at order 1/rriQ 
(loo = 2.0 ± 0.3 GeV, T = 0.8 ± 0.2 GeV). Since the order a s terms are large in the two-point 
function, to obtain a reliable result, it is imperative to include the order a s corrections into the 
analysis of the three-point functions as well. The evaluation of the arising two loop diagrams 
require somewhat sophisticated calculational techniques [23,20]. By now, all parameters that 
appear in the heavy quark expansion of meson form factors have been calculated to order 
a s /mQ in the framework of QCD sum rules* 

Isgur-Wise function: Several authors calculated the Isgur-Wise function [17,19-21] (see: Fig. 4) 
Of special interest is its slope at zero recoil, described by the p 2 parameter. It turns out, how- 
ever, that the sum rule prediction for this quantity is sensitive to the choice of the duality 
region, in particular to whether the continuum threshold ujq is allowed to depend on y [38]. 
Combining these uncertainties we conclude that QCD sum rules have a limited power in deter- 
mining p 2 , only constraining it to 0.6 < p 2 ^ 1.2 [17,19]. This upper limit can only be realized 
by allowing ujq to depend on y, and it slightly exceeds the Voloshin bound p 2 ^ 1. In fact, the 
experimental data seem to prefer p 2 « 1 [37]. 

1/rriQ corrections: For the four universal functions that appear at order 1/mg of the heavy 
quark expansion, the perturbative order a s corrections are even more significant compared to 
the leading order calculation [22], than in the case of the two-point function sum rule [23,24]. 
These corrections can change certain results by as much as a factor of two. The form factors 
X2 and X3) which parameterize the matrix elements of the O mag chromomagnetic operator in 
the effective Lagrangian, come out to be numerically negligible [23]. However, £3 is of order 
unity; it is proportional to the Isgur-Wise function to a very good approximation [24]. It 



★ This is not true for Xi(y)- However, this function does not affect any of the model-independent predictions 
of HQET. 
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Figure 4: Prediction for the Isgur-Wise function. The dashed lines indicate the Bjorken and 
Voloshin bounds on the slope at y = 1. The shaded region shows only the uncertainty related 
to variations of the sum rule parameters (from [13]). 



is convenient to present the result in terms of the renormalization group invariant ratio 

V(V) = |^4 = 0.6 ±0.2, (26) 

which already includes a conservative estimate of the systematic uncertainties as well. This 
is the result that we have already discussed in Section 2, which makes the extraction of |V^,| 
surprisingly accurate from the kinematically suppressed B — > D iv decay mode. 

In summary: QCD sum rules provide at present the theoretically most consistent framework 
to calculate the universal functions appearing in HQET, because: 

• The Ward identities of the effective theory, i.e. the zero recoil conditions are exactly 
fulfilled, independent of the sum rule parameters; 

• Due to the simple Feynman rules of HQET, two loop corrections can be included; 

• It allows for systematic renormalization group improvement; 

• There is a self-consistency check within the model on the assumption made. 

Some of the uncertainties can be estimated in a reliable manner. These are the ones related 
to (a) the truncation of the power series in the OPE (condensates); (b) higher orders in the 
perturbative expansion; (c) numerical values of the condensates. It is much harder to estimate 
the systematic uncertainties inherent in the method, related to the assumption of local duality. 
Therefore, it is crucial to determine carefully the sum rule window and check the stability, which 
(if found) makes us believe that duality indeed holds. Due to these systematic uncertainties, 
the accuracy of the QCD sum rule method cannot be arbitrarily improved (unlike the case of 



lattice calculations). Experience tells us that if all the above requirements are met, then the 
physical predictions of QCD sum rules are accurate to within 20-30%. 



4.3. Form factors 

The predictive power of the above results become more transparent if we translate them into 
predictions for the form factors measurable in B — > D^iV decays. Throughout the following 
discussion we use the form factor basis defined in Ref. [39]. In this notation, in the vanishing 
lepton mass limit, F\(q 2 ) is the single measurable form factor in B — > D Iv decays, while 
Ai(q 2 ), A2(q 2 ), and V{q 2 ) can be measured in B — > D^iv. In the infinite heavy quark mass 
limit [12] 

.2 n-l 



1 



Q 



(tub + m D *) 2 



W) = W)=W) = Fi(^). 



(27) 



The finite masses of the b and c quarks introduce symmetry breaking effects. The predictions 
for these form factors are plotted in Ref. [13]. It is useful to define ratios of the three form 
factors measurable in B — > D*£ v decays: 
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(28) 



Clearly, in the heavy quark symmetry limit R\ = R2 = 1. In the following table we present 
the predictions of our QCD sum rule analysis in HQET [23,24,13], together with the quark 
models of ISGW [40], BSW [39], KS [41], and QCD sum rules in "full" QCD [42] (not including 
radiative corrections): 



B — > D*£ v 


HQET + 
QCD SR 


ISGW 


BSW 


KS 


QCD SR 


#l(9max) 


1.35 


1.01 


0.91 


1.09 


1.31 


Ri(0) 


1.27 


1.27 


1.09 


1.00 


1.23 


^2(<7max) 


0.79 


0.91 


0.85 


1.09 


0.95 


ife(0) 


0.85 


1.14 


1.06 


1.00 


1.05 



Table 2: Predictions for the form factor ratios at the endpoins of the spectrum. 

The large values of R\ is an almost model-independent prediction of HQET, in the sense that 
both the QCD and the l/m c corrections are positive, and the 1/m^ terms could only cancel the 
previous two, if the QCD sum rule prediction for £3 failed with an order of magnitude. On the 
other hand, R2 depends more strongly on £3, and other models tend to give significantly higher 
values than ours. An experimental measurement of this quantity could distinguish between the 
above predictions. 



5. Summary and Open Problems 



In conclusion, we have seen how heavy quark symmetry reduces hadronic uncertainties by 
providing relations among heavy hadron form factors, and determining the absolute normal- 
ization of some of them at the kinematic limit point. Deviations from the predictions in this 
symmetry limit can be taken into account in the framework of a low energy effective theory. To 
calculate the universal functions appearing in HQET, QCD sum rules provide the theoretically 
most consistent framework. Of particular interest is the extraction of |V^,|: the theoretical 
uncertainty has been reduced to less than 5% in the decay B — > D*£ V, and it should be hardly 
larger in B — > D decays. Of course, the experimental measurement of B — > D iv near zero 
recoil is more difficult. However, it would be an independent measurement providing checks on 
both the theoretical and the experimental analysis. The QCD sum rule determination of the 
universal functions that appear at order 1/m.Q in the heavy quark expansion gives predictions 
for the symmetry breaking effects. A measurement of the form factor ratios in B — > D* 
decays would distinguish between these sum rule predictions and other models. 

In those cases when the final state does not contain a single heavy quark, HQS does not 
yield as restrictive relations as for heavy to heavy transitions. In particular, there is no absolute 
normalization of form factors. Still, different processes can be related to each other. This is 
the idea behind extracting \V u i\ from a comparison of B — > X iv and D — > X £&, where X is 
either n or p [43]. When X = ir, additional constraints are provided by chiral symmetry and 
soft pion theorems. The main problems are (a) that the contribution of the B* pole is large, 
and even dominant in the chiral limit; (b) the I/ttiq terms are expected to be significant. If 
these difficulties are overcome then the determination of \V u b\ from the experimental data may 
become much more reliable than that at present [44]. 

Here we only discussed exclusive semileptonic decays of heavy mesons. However, heavy 
quark symmetry also provides clues for accurate calculations of inclusive decays [36] . There are 
also exciting developments in progress to achieve a reliable, model-independent understanding 
of hadronic decays [45]. These theoretical developments will hopefully lead to a precise un- 
derstanding of heavy hadron physics, ultimately allowing to distinguish the signatures of new 
physics in heavy quark systems (if such signatures exist) from the ever decreasing uncertainties 
in the predictions of the Standard Model. 
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